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ABSTRACT 


A  high  resolution  20-level,  primitive  equation  model  is  initialized  using  sigma-0 
(a,)  data  acquired  during  a  Coastal  Transition  Zone  (CTZ)  cmise  from  6  to  12  July  1988 
near  Point  Arena,  Caliiomia  for  the  purpose  of  diagnosing  the  three-dimensional 
dynamically  balanced  flow  field  for  the  region.  The  major  feature  in  the  region  during 
the  cmise  period  was  a  strong  meandering  jet  which  flowed  equatorward  and  offshore, 
oriented  in  a  nonheast  to  southwest  direction.  A  density  diagnostic  version  of  the  model, 
in  which  the  field  is  held  constant  during  the  72  hour  integration  time,  is  compared 
with  a  robust  diagnostic  version  of  the  model,  in  which  the  field  is  allowed  to  adjust 
dynamically  within  the  constraints  of  a  damping  time  scale.  The  density  diagnostic 
version  of  the  model  produces  vertical  velocities  which  are  an  order  of  magnitude  larger 
than  vertical  velocities  calculated  from  bio-optical  data  collected  in  the  same  location  and 
at  the  same  time  as  the  CTZ  cmise.  The  robust  diagnostic  version  produced  vertical 
velocities  within  the  range  of  vertical  velocities  calculated  from  the  bio-optical  data.  The 
horizontal  pattern  of  vertical  velocity  from  the  robust  diagnostic  version  also  compared 
well  with  patterns  seen  in  trough/ridge  systems  in  the  Gulf  Stream.  In  particular,  there 
is  horizontal  convergence  and  downwelling  all  along  the  offshore  directed  part  of  the 
meandering  jet.  Vertical  cross-section  comparisons  between  model  velocity  data  and 
ADCP  velocity  data  collected  during  the  cmise  revealed  differences  which  may  be  caused 
by  factors  neglected  by  the  model  (e.g.  bottom  topography,  wind  forcing,  and  heat  fluxes). 
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1.  INTRODUCTION 


A.  BACKGROUND  -  THE  CALIFORNIA  CURRENT 

The  California  Current  is  the  southeastward  extension  of  the  North  Pacific  Drift 
which  flows  across  the  Pacific  Ocean  from  west  to  east  at  about  45°  N  (Tchemia,  1980). 
This  current  Cows  toward  the  equator,  and  is  considered  a  cold,  eastern  boundary  current. 
It.  along  with  three  other  currents,  comprises  the  California  Current  System.  The  other 
three  currents  are:  the  California  Undercurrent,  the  Davidson  Current  and  the  Southern 
California  Countercurrent  (Hickey,  1979).  These  currents  are  all  poleward  flowing,  the 
Davidson  Current  and  the  Southern  California  Countercurrent  being  surface  currents,  and 
the  California  Undercurrent  flowing  beneath  the  surface.  The  Davidson  Current  is  found 
north  of  Point  Conception  during  the  fall  and  winter.  The  Southern  California 
Countercurrent  exists  to  the  south  of  Point  Conception  in  the  California  Bight. 

The  California  Current  System  impacts  the  west  coast  states  climatically  and 
economically.  Upwelling  associated  with  it  brings  colder,  nutrient  rich  water  to  the 
surface,  supporting  marine  life  there  (Brink,  1983;  Huyer,  1983).  The  colder  water  at  the 
surface  also  cools  the  relatively  warm  air  flowing  over  it,  producing  mild  weather 
conditions  and  fog  along  the  coast  (Foster,  1989). 

Cold  filaments,  which  occur  in  the  California  Current  System  from  Vancouver 
Island  to  Point  Conception,  are  of  particular  interest.  These  filaments  consist  of  colil. 
salty  inshore  water  which  has  been  carried  as  far  as  300-400  km  offshore  from  the  coast 


(Strub  and  the  CTZ  Group,  1990).  Infrared  and  visible  color  satellite  images  reveal  these 
filaments  in  the  California  Current  System  as  well  as  in  other  eastern  boundary  currents 
(Traganza  et  al,  1983;  Abbott  and  Zion,  1985;  Flament  et  al.,  1985;  Kelly,  1985). 

The  Coastal  Transition  Zone  (CTZ)  experiment  examines  some  of  these  cold 
filaments  in  detail.  Its  purpose  is  to  study  the  physical,  chemical  and  biological  structure 
of  the  filaments  which  occuned  off  the  northern  California  coast  during  the  spring  and 
summer  of  1987  and  1988  (Strub  and  the  CTZ  Group,  1990).  The  goals  of  the  CTZ 
experiment  are;  understand  the  interaction  processes  between  the  cold  filament  water  and 
the  open  ocean  water;  study  the  effect  the  filaments  have  on  the  biology  and  metereology 
of  the  coastal  zone;  develop  maps  of  the  current,  hydrographic,  nutrient,  and  particle 
fields;  and  use  numerical  models  to  investigate  what  driving  forces,  dynamics,  topography 
and  stratiAcation  is  necessary  to  duplicate  the  observed  characterisitcs  of  these  filaments 
(CTZ  Group,  1988). 

B.  OBJECTIVE 

The  objective  of  this  thesis  research  is  to  diagnose  the  velocity  fields  obtained  from 
a  20-level  robust  primitive  equation  model  using  sigma-6  (Oe)  data  from  a  cmise  made 
during  the  1988  CTZ  field  experiment.  The  cruise  domain  is  from  approximately  37°  10' 
to  39°21'N  and  123°50'  to  126°40'W,  and  covers  about  one  meander  of  the  California 
Current. 
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The  same  domain  and  model,  lacking  the  robust  feature,  were  used  previously  to 
diagnose  the  velocity  fields  from  temperature  data  collected  during  the  same  experiment 
(de  Jesus,  1990). 
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II.  DATA  ACQUISITION  AND  PROCESSING 


A.  DATA  ACQUISITION 

Data  from  a  Naval  Postgraduate  School  cruise  made  from  6  to  12  July  1988  on  the 
R/V  POINT  SUR  that  was  part  of  the  CTZ  experiment  are  used  here.  This  cmise  is 
subsequently  referred  to  as  CTZRl. 

Acoustic  doppler  current  profder  (ADCP)  data  and  conductivity,  temperature  and 
depth  (CTD)  data  were  collected  during  the  cruise.  CTD  measurements  were  made  with 
an  NBIS  CTD  system.  Temperature  measurement  accuracy  and  salinity  measurement 
accuracy  were  about  ±0.01  C  and  ±0.003  psu  respectively.  Data  were  collected  to  a  depth 
of  500  m.  The  ADCP  made  by  RD  Instruments  operated  continuously  during  the  cruise. 
The  transducer  had  a  nominal  range  of  350  m  and  operated  at  150  kHz  with  a  4  m  pulse 
length.  Returning  echoes  were  sampled  selectively  in  4  m  bins  to  produce  a  vertically- 
averaged,  ship  relative  velocity.  Relative  velocity  profiles  were  made  every  three 
minutes.  Reliable  data  were  collected  to  depths  as  shallow  as  25  m  (Huyer  et  al.,  1990; 
de  Jesus,  1990). 

The  horizontal  sampling  grid  spacing  is  nominally  25  by  40  km  (Stanton  et  al., 
1990).  The  actual  station  grid  extends  approximately  120  nm  offshore  and  110  nm 
alongshore,  roughly  parallel  to  the  California  coast  and  southwest  of  Point  Arena  (Figures 
1  and  2).  The  alongshore  sections  are  approximately  20  nm  apart  and  have  intlivitlual 
stations  approximately  13  to  14  tun  apart  (de  Jesus,  1990).  Not  all  of  the  stations  in  the 
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CrZ  grid  are  included  in  the  model  domain  to  eliminate  coastal  effects  (Figure  2)  (de 
Jesus,  1990). 

B.  DATA  PROCESSING 

The  sigma-0  (o,)  data  were  derived  from  the  CTD  data  and  objectively  analyzed 
using  the  optimal  interpolation  method  described,  for  example,  by  Carter  and  Robinson 
(1987).  Tliis  method  and  some  of  the  following  procedures  were  similar  to  those  used 
by  de  Jesus  (1990)  for  temperature  data,  but  with  some  differences  as  described  below. 

The  spatial  mean  a,  profile  for  the  cruise  and,  the  objectively  analyzed  o,  at  each 
of  the  points  in  a  grid  with  4  km  horizontal  spacing  was  obtained  with  a  vertical  spacing 
corresponding  to  the  upper  7  levels  of  the  numerical  model  as  given  in  Table  1.  This 
mean  a,  profile  was  calculated  at  meter  depth  intervals  in  the  upper  500  ni,  corresponding 
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Figure  2.  Actual  Topography  in  the  Model  Domain:  The  box  shows  the  area  of  the 
model  domain.  Depths  are  in  meters. 

to  the  depths  at  which  observed  data  were  collected,  and  then  averaged  in  the  vertical  to 
every  10  m. 

The  mean  a*  profile  was  extended  to  4000  m  by  fitting  an  exponential  curve  to  the 
mean  profile  in  the  upper  500  m  using  a  centered  difference  scheme  around  450  m.  The 
O,  at  4000  m  was  left  to  be  detennined  by  the  exponential  curve,  increasing  the  degrees 


of  freedom  and  allowing  continuity  of  (jg  and  its  vertical  derivative  at  the  interface 
between  the  observations  and  the  extrapolation  (i.e.  at  450  m). 

This  extended  mean  o,  profile  was  used  to  calculate  the  first  three  dynamical 
vertical  modes  Ggilk),  GgiCk),  and  Ogjlk)  following  the  procedure  described  in  Rienecker 
et  al.  (1987)  and  outlined  below. 

The  initial  definition  made  was  the  perturbation 

Ggl/J,*)  =  <yg{ij,k)  -  Gg{^),  U) 

where  i,  j,  and  k  are  the  gridpoint  indices  for  the  x,  y,  and  z  directions,  respectively, 
G9'(i,jjc)  is  the  perturbation  part  of  Gg  at  the  gridpoint,  Ggli  j4c)  is  the  objectively  analyzed 
G,  at  the  gridpoint,  and  Gg(k)  is  the  mean  G,  at  level  k. 

To  extend  the  perturbation  below  500  m,  Ge'(i,j4c)  is  assumed  to  be  a  linear 
combination  of  the  first  three  dynamical  vertical  modes,  such  that 

=  a,{ijya^,{k)  +  a^[iJXy„{k)  +  aj{iJ)a„{k),  {2} 

where  a,(i,j)  is  the  weighting  given  to  Gei(k),  a2(i,j)  is  the  weighting  given  to  Ggjlk)  and 
aj(i,J)  is  the  weighting  given  to  GgjCk).  Solutions  for  a,,  a,,  and  a,  at  each  horizontal  grid 
point  ij  need  to  be  calculated. 

The  procedure  for  determining  a„  a^,  and  a,  followed  by  de  Jesus  (1990)  involved 
solving  a  system  of  four  equations  in  three  unknowns  (a,,  a^,  and  a,)  using  a  linear  least 
squares  method.  This  procedure,  referred  to  as  the  "mean  conserving”  fonnulation,  was 
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also  used  here.  In  addition,  two  new  procedures  were  used  which  gave  different  results 
for  the  vertical  profiles  below  500  meters.  They  were  considered  possibilities  for  bener 
representing  the  physical  condition  of  the  ocean  in  the  cruise  domain. 

The  "mean  conserving"  formulation  of  the  linear  least  squares  method  is  outlined 
briefly  as  are  the  two  new,  possibly  more  physical,  formulations.  For  the  sake  of 
comparison,  vertical  profiles  of  a*'  are  calculated  at  representative  (ij)  points  using  each 
of  the  three  methods.  Arguments  are  made  concerning  which  of  the  three  methods  are 
considered  the  best  for  representing  the  characteristics  of  the  ocean.  Finally,  an  exact 
method  of  solving  for  a,,  aj,  and  using  four  equations  and  four  unknowns  is  compared 
with  the  linear  least  squares  method,  and  arguments  are  made  for  using  the  linear  least 
squares  method  instead  of  the  exact  method.  These  issues  are  discussed  further  below. 

Following  the  "mean  conserving”  formulation  of  the  linear  least  squares  method, 
four  equations  are  derived  from  equation  (2): 

«Te(»,/^)>  =a,<o„(k)>  +aj<a„(*:)>  (^> 

«T„(A:)Oe(i\/,A:)>  =a,<a„(k)o„(k)>  +n2<a„(k)a„(k)>  *a^<0,^{k)0„{k)>, 
«J„{kWf{iJ,k)>  =a,«y„(k)0„(A:)>  ■>-a^«y^{k}o„{k)>  +a,<a„(A:)ag3(k)>,  (5} 

«J„{kW,{iJ,k)>  =a,«y,3(k)0„(k)>  *a^«J„{k)o„{k)>  *a^«3„{kXj^^{k)>.  (<») 

The  bracketed  (<  >)  teniis  are  a  vertical  average,  defined  as 

7  7 

<Q>  s  (2  A  z{k)  /  Y,  A 

t  I 
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where  Q  is  any  expression  inside  the  brackets,  Az{k)  is  the  vertical  thickness  between 
adjacent  k  levels  in  the  model,  and  the  limits  1  and  7  with  the  summation  sign  correspond 
to  levels  I  and  7  in  the  model(Table  1).  Therefore,  <Q>  is  the  vertical  average  of  Q(k) 
over  the  top  seven  levels  of  the  model,  which  are  encompassed  by  the  analyzed 
observations  (i.e.  the  upper  500  meters).  Equation  |3)  is  the  statement  of  "conservation 
of  the  mean"  in  this  formulation. 

This  system  of  four  equations  in  three  unknowns  (a„  a^,  a,)  is  solved  by  the  linear 
least  squares  method.  This  means  that  none  of  the  four  equations  are  satisfied  exactly, 
but  that  there  is  a  residual  error  in  each  one.  After  a„  aj,  and  a,  are  found  for  each  i,j 
point  in  the  domain  of  the  model,  equation  |2)  is  used  to  find  o,'(i,j4c)  for  k  =  8  to  k  = 
20  (i.e.  from  500  meters  to  4000  meters)  (Table  I).  From  Equation  { 1 ),  a9(i,jjk)  is  found 
for  all  20  levels  in  the  model. 

The  first  of  the  two  new  formulations  alluded  to  above  represents  an  energy-related 
approach  to  arriving  at  solutions  for  a,,  aj,  and  a^.  It  involves  the  replacement  of  equation 
|3)  by 

<oi(f  j»oi(i  J,*)>  »a,<oi(i\/,*)0„(*)>  +a,<oi(i\/, *)«„(*)> 

The  on  the  left  hand  side  of  the  equation  is  an  energy  term.  Since  this  equation  puts 
a  constraint  on  the  energy  of  the  in  the  upper  ocean  where  the  observations  occur,  this 
is  called  the  "energy  conserving"  formulation  of  the  problem. 


9 


The  second  of  the  new  formulations  puts  a  constraint  on  the  continuity  of  at  the 
interface  between  the  observations  and  the  region  of  extrapolation  below.  In  this  case, 
the  interface  occurs  at  level  7  in  the  model.  Therefore,  equation  { 3 }  is  replaced  by 

=a,<Te,(7)  (9} 

Because  of  the  constraint  at  the  interface,  this  is  called  the  "interface  matching" 
formulation.  Since  both  of  the  new  formulations  are  also  solved  using  a  linear  least 
squares  method,  there  is  a  residual  error  in  each  equation  as  there  is  with  the  first 
formulation. 

When  comparing  the  results  of  vertical  profiles  using  these  three  formulations, 
differences  among  them  can  readily  be  seen.  In  the  CTZRl  cmise  domain,  two  points 
where  o*'  is  large  in  the  upper  levels  of  the  ocean  are  used  to  illustrate  these  differences 
(Figures  4  and  5). 

The  (i,j)  gridpoint  (27,36),  which  is  qjproximately  mjq)  coordinate  (x,y)  point 
(1 10,144)  in  Figure  3,  is  located  in  a  region  of  low  density  water,  with  <  24.4  kg/m^ 
at  65  m.  Figure  4  shows  the  vertical  profile  of  a/  at  this  point  for  the  three  different 
formulations  of  the  linear  least  squares  solution  to  a„  aj,  and  a,  in  equation  1 2 ) .  Clearly, 
the  observed  profiles  in  the  upper  500  meters  are  the  same  in  each  case.  Also,  below 
about  2000  meters  the  proHles  are  the  same.  In  the  deep  ocean,  the  water  characteristics 
are  such  that  the  dynamical  modes  are  all  small  in  amplitude  so  that  changing  the 
formulation  of  the  solution  does  not  affect  the  profile.  Between  these  two  levels,  the 
three  foimulations  produce  three  different  profiles.  The  a,'  profile  for  the  "interface 
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Figure  3.  Plot  of  o,  at  65  Meters  Depth(Model  Level  3),  Cruise  CTZRl,  06-12 
July  1988:  The  two  points  of  interest  are  marked. 

matching"  formulation  has  the  smallest  amplitude  of  all.  The  constraint  on  having 
continuity  at  the  interface  between  the  observations  and  the  extrapolation  is  the  cause  of 


this  reduced  amplitude.  The  "energy  conserving"  fonnulation  is  also  smaller  ui  amplitude 
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DEPTH  IN  METERS 


disturbrnce:  sigmr  0  in  kg/m' 


Figure  4.  Vertical  Profiles  of  a,'  at  Map  Coordinate  Point  (110,144): 

The  symbols(located  at  the  20  model  levels)  and  the 
lines  for  each  fonnulation  are  given  by  the  legend. 
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than  the  "mean  conserving"  formulation  of  the  solution.  Comparing  the  profiles  in  Figure 
4  with  the  three  vertical  dynamical  modes  shown  in  Figure  6  suggests  that  mode  2  and 
mode  3  are  the  main  contributors  to  the  stmcture  of  the  vertical  profile  of  at  the  level 
between  500  and  2000  meters. 

The  (i,j)  gridpoint  (21,8),  which  roughly  corresponds  to  map  coordinate  (x,y)  point 
(85.32).  is  in  an  area  of  high  density  water  at  65  m  as  shown  in  Figure  3.  The  vertical 
profiles  of  the  three  formulations  at  this  point  are  not  as  distinctly  different  as  those  at 
point  (1 10,144)  as  seen  when  comparing  Figure  5  with  Figure  4.  At  point  (85,32),  the 
"interface  matching"  formulation  and  the  "energy  conserving"  formulation  have  a  smaller 
amplitude  than  the  "mean  conserving"  formulation  in  the  level  between  500  meters  and 
2000  meters.  This  illustrates  the  differences  among  the  points  in  the  cruise  domain 
resulting  from  the  3  different  extrapolation  formulations.  On  the  other  hand,  the 
siiTularities  of  the  profiles  above  500  meters  and  below  2000  meters  can  be  explained 
using  the  arguments  made  previously  for  point  (110,144). 

Of  the  three  formulations  described  above,  which  is  the  "best"  to  use  for  input  into 
the  model?  The  "mean  conserving"  formulation  suffers  because  it  places  no  constraint 
on  the  amplitude  of  The  best  example  of  this  is  a  sine  wave.  The  mean  of  a  sine 
wave  is  the  same  no  matter  what  its  ampliturte.  The  "energy  conserving"  formulation 
may  alias  high  wave  number  energy  observed  in  the  upper  500  m  (which  actually  belong 
to  higher  dynamical  modes)  onto  the  lower  wave  number  modes,  which  is  due  to  the 
energy  dependence  caused  by  the  squared  term  on  the  left-hand  side  of  equation  ( H 1 
This  would  tend  to  increase  the  amplitude  of  a,'  compared  to  a  method  that  did  not  alias 
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the  high  wave  number  variability  into  the  solution.  Therefore,  the  "interface  matching" 
formulation  seems  to  be  the  "best"  to  use  for  the  model.  As  stated  above,  it  places  a 
constraint  on  the  <5^  at  the  interface  between  the  observed  data  and  the  extrapolated  data, 
and  does  not  alias  high  wave  number  variability  onto  the  three  low  wave  number  modes. 
The  problems  associated  with  the  other  two  methods  are,  therefore,  eliminated.  However, 
before  a  definite  conclusion  can  be  made  about  the  relative  ability  of  these  formulations 
to  depict  sub-500  m  variability,  they  should  be  tested  with  full  water  column  observed 
data.  Because  of  the  arguments  given  above,  the  "interface  matching"  formulation  is  used 
in  this  study.  Maps  of  a,,  aj,  and  83  using  the  "interface  matching"  formulation  are  shown 
in  Figures  7,  8,  and  9.  These  modal  amplitudes  have  large  spatial  coherence,  and  their 
patterns  closely  correspond  to  the  o,  field  in  Figure  3. 

The  problem  of  solving  for  a,,  a,,  and  a,  in  equation  |2)  could  have  been  done 
exactly  with  four  equations  and  four  unknowns  using  variational  techniques  found  in 
Hildebrand  (1965).  In  Hildebrand  (1965;  see  pp.  119-121),  it  is  shown  that  the  problem 
of  determining  a,,  and  a,  that  minimizes  the  root  mean  square  difference  between  the 
observed  o,'  and  the  <j,'  given  by  (2)  in  the  upper  500  m  subject  to  the  constraint  (9), 
can  be  expressed  as  an  unconstrained  minimization  problem  which  leads  to  the  following 
four  equations 
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Figure  7.  Plot  of  A,  Values,  Cruise  CTZRl:  Obtained  from  "interface  matching" 
formulation  described  in  text. 
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Figure  8.  Plot  of  Aj  Values,  Cruise  CTZRl;  Obtained  from  "interface  matching" 
formulation  described  in  text. 


Figure  9.  Plot  of  A,  Values,  Cruise  CTZRl:  Obtained  from  "interface  matching" 
formulation  described  in  text. 
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{10} 


=dr,0„(7)  +«jO,j(7)  +flja„(7), 

<0,^{k)o't{ij,k)>  -a,«J„(k)0^^{k)>  +dj<0„(i)<J„(t)>  +a,<0„(fc)a,j(*)>  - Aog,(7),  {11} 

\ 

«Jgj(l:)Og(/V,A:)  =a,<0gj(*)0„(*)>  +aj<a,j(*)0,j(fc)>  *a^«J„(k)C!„{k)>  -  —0,^(7),  {12} 

A, 

<0„{k)o',iiJ,k)>  =a,<0„(*)0„(^)>  +<Jj<a,j(*)o,j(*)>  +03«Jg3(*)<J„(*)> -—a„0),  {13} 

where  X  is  the  (unknown)  Lagrange  multiplier.  Equation  {10)  is  the  same  as  equation 
(9).  Equations  (11),  (12),  and  (13)  are  the  same  as  equations  (4),  (S),  and  (6), 
respectively,  with  the  exception  of  the  final  term  on  the  right  hand  side  of  the  equations. 
This  flnal  term,  the  term  which  includes  the  Lagrange  multiplier,  is  a  residual  term. 
Notice  that  this  exact  method  allows  the  constraint,  in  this  case  (10),  to  be  satisfied 
exactly.  All  the  residual  error  is  put  into  the  other  3  equations.  As  pointed  out  above, 
the  least  squares  method  distributes  the  error  among  all  4  equations,  none  of  which  are 
satisfied  exactly.  These  residuals  represent  the  extent  to  which  a*'  as  represented  by  the 
first  three  vertical  dynamical  modes  does  not  match  the  actual  ae'  in  the  region  of  the 
ocean  where  actual  observations  have  been  made,  in  this  case  the  upper  500  m.  The 
residual  term,  therefore,  represents  a  measure  of  the  high  wave  number  variability  above 
500  m  that  simply  cannot  be  represented  by  a  sum  of  the  first  3  vertical  dynamical 
modes.  Since  o,'  may  have  high  wave  number  variability  at  all  levels  at  and  above  500 
m,  it  does  not  seem  physically  realistic  to  force  an  exact  contraint  at  level  7  of  the  model 
as  in  equation  (10). 
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Thus,  die  least  squares  method,  with  4  equations  and  3  unknowns,  was  chosen  for 
this  study.  In  this  least  squares  method,  none  of  the  equations  was  satisfied  exactly,  but 
the  total  residual,  spread  over  all  4  equations,  may  be  smaller  than  it  would  be  using  the 
exact  method.  This  is  also  a  subject  worthy  of  much  further  study. 

C.  THE  MODEL 

The  numerical  model  is  a  20-level,  primitive  equation  model  of  a  baroclinic  ocean 
on  an  f-plane,  with  hydrostatic,  Boussinesq,  and  rigid  lid  approximations,  and  a  flat 
bottom  (Haney,  1985).  The  space-staggered  B-scheme  is  used  for  horizontal  finite 
differencing,  and  a  sigma  coordinate  system  defines  the  vertical  dimension  (Arakawa  and 
Lamb,  1977).  The  boundary  conditions  at  the  sea  surface  consist  of  no  wind  stress  and 
no  buoyancy  flux.  The  lateral  boundaries  of  the  model  domain  are  treated  as  open 
boundaries.  A  sponge  layer  of  5  giidpoints  based  on  the  discussion  of  Lorenzetti  and 
Wang  (1986)  is  used  at  the  lateral  boundaries  around  the  entire  model  domain  to 
eliminate  the  effects  of  gravity  wave  noise  diere.  The  coefficients  of  the  sponge  at  the 
gridpoints,  starting  from  die  boundary  inward,  are  .20,  .13,  .07,  .03,  and  .008, 
respectively.  The  sponge  affects  only  the  velocity  fields  in  the  model.  The  total  depth 
of  the  ocean  is  4000  meters.  Model  levels  are  listed  in  Table  1 .  The  u  and  v  components 
of  velocity,  and  o,  are  calculated  at  the  model  levels,  but  the  w  component  of  velocity 
is  calculated  at  the  midpoints  between  the  model  levels  as  a  vertically  staggered 
arrangement.  The  values  of  constants  used  in  the  model  are  listed  in  Table  2. 
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The  CTZRl  cruise  domain  included  76  objectively  analyzed  gridpoints  in  the  east- 
west  direction  by  60  gridpoints  in  the  north-south  direction.  The  model  domain  included 
only  64  of  the  gridpoints  in  the  east-west  direction  by  60  gridpoints  in  the  north-south 
direction,  with  the  12  eliminated  gridpoints  in  the  east-west  direction  being  in  the 
eastern  part  of  the  cruise  domain  along  the  California  coast.  This  makes  the  model 
domain  more  uniform,  and  reduces  any  shallow  water  effects  which  may  occur. 

Since  the  model  grid  points  are  much  closer  together  than  the  observed  data  points, 
the  objectively  analyzed  a,  field  was  somewhat  noisy  on  the  fine  4  km  model  grid.  Since 
this  noise  is  not  considered  to  be  physical,  the  following  horizontal  smoother  was  used: 

where  a«'(i,j4£),  the  perturbation  at  each  domain  gridpoint, 

O,*  =  cr,'(i+l,j4c)  +  Oe'(i-l  j4c)  +  <T,'(i,j+lJk)  +  0,'(i,j-14c),  and 

a,*  =  Oe^(i-Kl,j+14c)  +  o,'(i-l  j-^14c)  +  o/(i-l  j-ljc)  +  a0'(i+l,j-14c).  This  smoother  was 
applied  twice.  Figure  3  is  an  example  of  a  smoothed  a,  field.  The  corresponding 
unsmoothed  field  is  shown  in  Figure  10. 

Two  methods  were  used  to  diagnose  the  currents  from  a  given  o,  field.  In  both 
methods,  the  c,  field  and  the  geostrophic  currents  (having  zero  vertical  average  over  the 
full  model  depth)  computed  from  the  hydrostatic  pressure  field  were  used  to  initialize  the 
20  level  model.  The  model  was  then  integrated  forward  in  time  (for  several  days)  until 
the  currents  and  density  (a,)  fields  came  into  a  new  "balance."  In  one  method,  the 
density  diagnostic  method,  the  o,  field  is  held  constant  throughout  the  model  run.  thus 
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Table  1.  K-LEVELS  AND  CORRESPONDING  DEPTHS  IN  METERS  USED  IN 
MODEL. 


k  level 

z(k)  (meters) 

k  level 

z(k)  (meters) 

1 

5 

11 

1105 

2 

25 

12 

1325 

3 

65 

13 

1565 

4 

125 

14 

1825 

5 

205 

15 

2105 

6 

305 

16 

2405 

7 

425 

17 

2725 

8 

565 

18 

3065 

9 

725 

19 

3425 

10 

905 

20 

3805 

requiring  the  velocity  fields  to  adjust  to  it.  In  the  second  method,  a  robust  diagnostic  was 
added,  in  which  the  <!•  Held  was  allowed  to  adjust,  but  was  damped  back  to  the  initial 
features  on  a  time  scale  which  was  short  compared  to  the  advective  time  scale,  but  long 
compared  to  the  geostrophic  adjustment  time  scale  (Saimiento  and  Bryan,  1982).  This 
was  done  to  hold  the  larger  scale  features  constant,  therefore,  maintaining  the  diagnostic 
feature  of  the  model,  but  allowing  the  smaller  scale  features  to  adjust.  Much  of  the  small 
scale  noise  arises  from  the  multi-level  initialization  from  the  objectively  analyzed,  sparse 
observations.  The  scaling  arguments  used  to  justify  the  choice  of  the  damping  time  scale 
are  presented  in  the  results  section  below. 
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in.  MODEL  RESULTS  AND  DISCUSSION 


The  primitive  equation  model  was  initialized  with  the  a,  fields  at  the  levels 
described  above,  and  with  geostrophically  balanced  currents  having  a  zero  vertical 
average.  The  model  was  then  integrated  forward  in  time  for  72  hours.  In  the  density 
diagnostic  version  of  the  model,  the  field  was  held  fixed.  During  this  integration,  the 
currents  underwent  adjustments  to  a  new  state  of  balance  in  which  non-linearity  and  eddy 
diffusion  effects  were  not  negligible  because  the  currents  were  no  longer  in  exact 
geostrophic  balance.  Time  series  of  the  velocity  components  are  used  to  show  the  process 
by  which  this  balance  occurs.  The  resulting  diagnostic  currents  are  discussed  below. 

In  addition,  the  robust  diagnostic  version  of  the  model,  as  briefly  described  above, 
was  integrated  forward  in  time  to  72  hours.  The  effects  of  the  robust  diagnostic  version 
versus  the  density  diagnostic  version  are  discussed.  Special  attention  is  given  to  the 
vertical  component  of  the  velocity  field  and  its  relationship  to  vorticity  advection  for  the 
robust  diagnostic  version  of  the  model.  Finally,  a  con^aiison  is  made  between  the 
currents  as  computed  by  the  robust  diagnostic  version  of  the  model  and  observed  currents 
measured  witfi  an  ADCP  during  CIZRl. 

To  simplify  the  presentation  and  discussion  of  the  results,  only  horizontal  maps  and 
time  series  at  65  m  depth  are  used,  which  is  a  model  level  representing  surface  currents 
in  the  thermocline.  Discussion  of  change  in  the  vertical  is  limited  to  the  comparison 
between  ADCP  results  and  model  results. 
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Table  2.  CONSTANTS  USED  IN  MODEL 


Constant 

Name 

Value 

n 

earth  rotation  rate 

2rcday' 

DTAU 

timestep 

450  s 

dx 

meridional  grid  spacing 

4  X  10’  cm 

dy 

zonal  grid  spacing 

4  X  lO’  cm 

D 

total  ocean  depth 

4  X  lO*  cm 

LAT 

reference  latitude 

38.25°  N 

f 

Coriolis  parameter 

9.00  X  10  ’  sec  ' 

g 

acceleration  of  gravity 

980  cm  sec'^ 

Am 

biharmonic  momentum  dif¬ 
fusion  coefficient 

4  X  10'^  cm^  sec  ' 

A„ 

biharmonic  heat  diffusion 
coefficient 

4  X  10'^  cm^  sec  ' 

Km 

vertical  eddy  viscosity 

0.5  cm*  sec ' 

K„ 

vertical  eddy  conductivity 

0.5  cm*  sec  ' 

a 

thermal  expansion  coeffi¬ 
cient 

2.877  X  10-*  C 

A.  TIME  SERIES  OF  VELOCITY  FIELDS  FOR  DENSITY  DIAGNOSTIC 
Was  a  new  state  of  balance  achieved  by  integrating  the  density  diagnostic  version 
of  the  model  ahead  72  hours?  Time  series  of  the  u,  v,  and  w-components  of  velocity  at 
(ij)  gridpoint  (44,25)  (i.e.  map  coordinate  point  (172,96)  in  Figure  11)  and  at  a  depth  of 
65  meters  (90  meters  for  w  because  of  the  vertically  staggered  grid  described  in  Chapter 
II)  show  the  adjustments  in  the  velocity  at  these  locations  (Figures  12  and  13).  Note  the 
large  oscillations  which  occur  during  the  early  integration  times  in  the  model.  These 
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Figure  10.  Unsmoothed  plot  of  a,. 

oscillations  ate  damped  as  the  new  state  of  balance  is  approached  at  around  48  hours. 
Tlie  stnaU  oscillations  continue  to  occur  at  72  hours,  especially  in  the  u  and  v  component 
velocities  (Figure  13).  The  oscillations  are  reninants  of  the  initial  adjustments  in  the 
velocities,  and  could  possibly  continue  for  several  more  days. 
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This  process  can  also  be  seen  in  the  time  series  of  the  u,  v,  and  w-con^nents  of 
velocity  at  (i,j)  gridpoint  (41,24),  or  (x,y)  map  coordinate  point  (160,92)  in  Figure  11. 
also  at  a  depth  of  65  meters  (90  meters  for  w)(Figures  14  and  15).  At  this  point,  the 
damped  oscillations  beyond  48  hours  in  the  u  and  v  velocity  components  are  almost 
sinusoidal. 

Examination  of  the  time  series  of  the  component  velocities  at  these  two  points  show 
that  an  exact  balance  has  not  been  reached  in  the  density  diagnostic  version  of  the  model 
at  72  hours.  Small  oscillations  are  still  ocurring,  although  it  appears  that  the  currents  are 
quite  close  to  their  final  equilibrium  values  at  72  hours.  The  similar  evolution  of  the 
velocity  components  with  time  at  the  two  points  indicates  that  the  entire  velocity  field 
probably  evolves  in  a  similar  way.  Integrating  the  model  beyond  72  hours  to  determine 
when  the  oscillations  in  the  velocity  components  are  fully  damped  was  deemed 
unnecessary. 

B.  SIGMA-e  AND  VELOCITY  FIELDS  FOR  DENSITY  DIAGNOSTIC 

The  <j,  field  at  time  0  for  a  depth  of  65  meters  is  shown  in  Figure  16,  and  the 
resulting  hydrostsatic  pressure  field  in  Figure  11.  By  definition,  these  fields  do  not 
change  with  time  in  the  density  diagnostic  model. 

The  initial  (geostrophic)  u  and  v  components  of  the  velocity  field  at  65  meters  are 
shown  in  Figures  17  and  18.  The  resulting  vector  velocity  field  is  shown  in  Figure  19. 
Note  the  strong  south  to  southwestward  jet  flowing  from  the  northeastern  comer  to  the 
southwestern  comer  of  the  model  domain.  Velocities  in  this  jet  reach  as  high  as  8(1 
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Cruise  CTZRI,  06-12  July  1988 

P  hour  7Z  depth  65  m 


Distonce  (km) 


Figure  11.  Pressure  Field  of  the  Density  Diagnostic  at  65  Meters  and  72  Hours: 

Points  referred  to  are  marked.  Contour  interval  2.5  cm.  Dashed 
lines  are  negative  values. 

cm/sec.  Velocities  as  high  as  75  cm/sec  were  observed  at  50  meters  by  Huyer  et  al. 
(1990)  during  CTZRI.  The  velocities  shown  in  the  model  domain  at  time  0  are  simply 
the  geostrophic  velocities  with  wliich  the  model  is  initialized.  The  initial  w  velocity 
component  is  zero  by  definition. 
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TIME  SERIES  OP  W 


TIME(HRS) 


Figure  12.  Time  Series  of  W  Velocity  Component  at  Map  Point  (172,96)  and 
Depth  of  90  Meters  for  Density  Diagnostic. 


After  integrating  the  model  ahead  72  hours,  the  velocities  are  no  longer  geostrophic. 
The  u  component  velocity  field  has  intensified,  especially  in  the  south  centra!  pait  of  the 
model  domain  where  the  jet  turns  wcstward(Figuic  20).  Tlie  v  component  velocity  field 
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TIME  SERIES  OF  U  AND  V 


TIME(HRS) 


Figure  13.  Time  Series  of  U  and  V  Velocity  Components  at  Map  Point  (172,96) 
and  Depth  of  65  Meters  for  Density  Diagnostic. 


has  also  adjusted(Figute  21).  This  intensiflcation  and  adjustment  can  also  be  seen  in  the 
time  series  shown  in  Figure  13.  Note  the  inciea.sed  u  velocity  shown  at  72  hours  as 
compared  to  0  hours.  Also  note  how  the  v  velocity  has  changed. 
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TIME  SERIES  OF  W 


Figure  14.  Time  Series  of  W  Velocity  Component  at  Map  Point  (160,92)  and 
Depth  of  90  Meters  for  Density  Diagnostic. 

The  resultant  horizontal  velocity  vectors  for  72  hours  are  shown  in  Figure  22.  A 
comparison  of  a  few  of  the  vectors  shown  here  with  the  corresponding  initial  geostrophic 
velocities  in  Figure  19  reveals  significant  direction  shifts  as  well  as  speed  changes.  The 
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TIME  SERIES  OF  U  RND  V 
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Figure  15.  Time  Series  of  U  and  V  Velocity  Components  at  Map  Point  (160,92) 
and  Depth  of  65  Meters  for  Density  Diagnostic. 


direction  shift  varies  from  point  to  point  as  illustrated  by  the  time  series  in  Figures  13  and 
15.  From  the  u  and  v  components  of  velocity  in  Figure  13,  it  can  be  ii\ferred  that  tlie 
current  has  shifted  westward  by  72  hours.  The  direction  of  the  cuiient  at  the  poutt  in 
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CRUISE  CTZRI,  06-12  JULY  1988 

SIGMA  0  HOUR  0.  DEPTH  65  M 


Figure  16.  Sigma-S  Field  at  Time  0  and  Depth  of  65  Meters:  Contour  interval  .2 


kg/m’ 

Figure  15  seems  to  have  shifted  to  the  east. 

The  maximum  velocities  shown  in  the  jet  at  72  hours  arc  over  100  cm/sec.  This 
is  higher  than  the  maximum  observed  values  noted  above. 

The  w  component  of  velocity  at  72  hours  and  90  meters  is  shown  in  Figure  23. 
Maximum  upwelling  is  around  .8  cm/sec,  and  maximum  downwcUuig  is  around  I O 
cm/sec.  At  the  points  shown  in  the  above  time  series,  upwelling  is  occurring  (Figures  12 
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Figure  17.  U  Component  Velocity  Field  at  Time  0  and  Depth  of  65  meters: 

Contour  interval  20  cm/s.  Dashed  lines  are  negative  values. 

and  14).  At  map  point  (172,96)  upwelling  velocity  is  around  .25  cm/sec  (Figure  12),  and 
at  map  point  (160,92)  upwelling  velocity  is  around  .6  cm/sec  (Figure  14),  at  72  hours. 
Tliesc  velocities  convert  to  200  to  700  meters  per  day  which  are  much  higher  than  tlie  6 
to  25  meters  per  day  observed  by  Washburn  et  al.  (1990)  and  Kadko  et  al.  (1990)  in  the 
CTZ  region. 


Figure  18.  V  Component  Velocity  Field  for  Time  0  and  Depth  of  65  Meters: 

Contour  interval  20  cin/s.  Dashed  lines  arc  negative  values. 

Therefore,  the  density  diagnostic  version  of  tlie  model  does  not  give  horizontal  or 
vertical  velocities  which  agree  closely  with  those  observed.  One  possible  reason  for  the 
rather  poor  results  for  the  density  diagnostic  method  is  that  the  currents  must  adjust 
entirely  to  tire  given  density  field.  However,  in  the  real  ocean  there  is  a  continual  mutual 
adjustment  between  tire  currents  and  the  density.  Futhermore,  on  scale.s  of  motion  that 
are  smaller  than  the  internal  Rossby  radius  (about  30  km  ui  the  CTZ  region)  the  dcirsiiy 


Figure  19.  Resultant  Horizontal  Velocity  Vectors  at  Time  0  and  Depth  of  65 
Meters:  Vectors  are  plotted  at  every  4th  gridpoint. 

field  would  actually  adjust  to  the  currents.  Thus,  the  density  diagnostic  method  is  not 
really  valid  for  scales  of  motion  of  the  order  of,  or  smaller  than,  about  30  km.  In  order 
to  facilitate  a  mutual  geostrophic  adjustment  between  the  currents  and  the  density  field, 
the  robust  diagnostic  method  is  introduced  in  the  next  section. 
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Figure  20.  U  Component  Velocity  at  Time  72  Hours  and  Depth  of  65  Meters  for 
Density  Diagnostic:  Contour  interval  20  cm/s.  Dashed  lines  are 
negative  values. 

C.  THE  ROBUST  DIAGNOSTIC  AND  ITS  EFFECTS 

The  robust  diagnostic  method,  as  described  in  Chapter  II,  allows  the  currents  and 
the  a,  field  to  mutually  adjust  to  each  other.  The  field  is  damped  back  to  its  initial 
value  on  a  time  scale  that  is  fast  compared  to  advection  but  slow  compared  to  the  time 
needed  for  geostropliic  adjustment.  As  stated  above,  this  type  of  damping  allows  (he 
smaller  scale  features  that  have  short  time  scales  to  mutually  adjust,  but  it  holds  the  larger 
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Figure  22.  Resultant  Horizontal  Velocity  Vectors  at  Time  72  Hours  and  Depth  of 
65  Meters  for  Density  Diagnostic:  Vectors  are  plotted  at  every  4th 
gridpoint. 

where  L  and  U  are  the  length  scale  and  the  velocity  scale,  respectively.  In  this  case,  L 
was  chosen  as  50  km  and  U  was  chosen  as  50  cm/s  which  yields  a  t,  of  approximately 
10’  seconds  or  1  day. 

The  geostrophic  adjustment  time  scale  is: 


CRUISE  CTZRl,  06-12  JULY  1988 

W  HOUR  72.  DEPTH  90  M 


DISTANCE  IN  KM 


Figure  23.  W  Component  Velocity  at  Time  72  Hours  and  Depth  of  90  Meters  for 
Density  Diagnostic:  Contour  interval  .2  cm/s.  Dashed  lines  are 
negative  values. 


wliere  f  is  the  coriolis  parameter.  For  this  case,  t,  is  approximately  10*  seconds  or  1/10 
of  a  day. 

The  damping  time  scale,  therefore,  should  be  much  longer  than  1/10  of  a  day  and 
much  shorter  than  1  day  to  produce  the  desired  effect.  Clearly,  it  is  not  possible  to 
satisfy  both  of  these  constraints.  For  this  study,  1  day  (t.)  was  chosen  as  the  damping 
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time  scale  to  maximize  the  time  over  which  geostrophic  adjustment  can  occur.  Since  the 
damping  time  is  the  same  as  the  advective  time,  and  not  smaller,  some  advective  effects 
will  be  unavoidable. 

Time  series  of  a,,  and  the  u,  v,  and  w  components  of  velocity  are  used  to  illustrate 
the  adjustment  process  with  the  robust  diagnostic  (Figures  24,  25,  26,  27,  28,  and  29). 
The  time  series  of  o,  for  map  points  (172,96)  and  (160,92)  and  a  depth  of  65  meters 
shown  in  Figures  24  and  27  indicate  a  downward  trend  with  time  of  a®  occurring  at  both 
points  for  the  robust  diagnostic  method. 

The  time  series  of  the  component  velocities  show  an  initial  adjustment  as  in  the 
density  diagnostic  case  followed  by  a  series  of  long  oscillations  (Figtues  25,  26,  28,  and 
29).  These  long  oscillations  seem  to  begin  at  about  30  hours  for  the  u  and  v  component 
velocities,  and  somewhat  later  for  the  w  velocity  component.  The  oscillation  is 
particularly  flat  for  the  v  velocity  component  at  map  point  (160,92).  The  long  oscillations 
at  both  points  are  probably  caused  by  advective  effects  at  the  longer  integration  times,  due 
to  using  1  day  (t.)  as  the  danping  time  scale.  In  other  words,  the  longer  scale  synoptic 
features  are  not  being  held  perfectly  constant.  At  the  model  grid  point  shown,  there  is 
sinking  motion  (w  <  0  in  Figure  25)  on  the  longer  time  scale.  The  resulting  vertical 
advection  of  the  mean  stratification  causes  the  density  (a«)  to  decrease  with  time  at  65 
m  (Figure  24).  The  model  is  losing  its  strictly  diagnostic  character.  This  trade-off  is, 
however,  desirable  in  terms  of  other  features  of  the  robust  diagnostic  as  discussed  below. 

Comparisons  between  the  u  and  v  velocity  component  time  series  for  the  density 
diagnostic  (Figures  13  and  15)  and  those  of  the  robust  diagnostic  (Figures  26  and  29) 
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TIME  SERIES  OF  SIGMR  0 


indicate  that  the  adjustment  in  horizontal  velocity  occurs  faster  for  the  robust  diagnostic 
case.  The  short  oscillations  which  are  characteristic  of  the  density  diagnostic  method  aie 
not  present  in  the  robust  diagnostic  case.  Whereas  the  density  diagnostic  has  not 


Figure  25.  Time  Series  of  W  Velocity  Component  at  Map  Point  (172,96)  and 
Depth  of  90  Meters  for  Robust  Diagnostic. 

completely  adjusted  by  72  hours,  the  robust  diagnostic  has,  discounting  the  long  wave 
oscillations  which  are  believed  to  be  advective  rather  than  an  initial  adjustment. 
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TIME  SERIES  OF  U  RND  V 
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Figure  26.  Time  Series  of  U  and  V  Velocity  Components  at  Map  Point  (172,96) 
and  Depth  of  65  Meters  for  Robust  Diagnostic. 

Coinparisons  between  the  density  diagnostic  and  robust  diagnostic  time  series  of  the 
w  component  velocity  may  be  misleading  (Figures  12,  14,  25.  and  28).  Tlic  lobu.si 
diagnostic  w  component  velocity  time  series  seem  to  be  oscillating  much  more  than  their 
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TIME  SERIES  OF  SIGMR  0 
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Figure  27.  Time  Series  of  a,  at  Map  Point  (160,92)  and  Depth  of  65  Meters  for 
Robust  Diagnostic. 

density  diagnostic  counterparts.  Tlie  robust  diagnostic  w  component  velocities  are, 
however,  an  order  of  magnitude  smaller  than  those  of  the  density  diagnostic.  Therefore, 
what  may  seem  like  larger  oscillations  are  not.  Tlie  order  of  magnitude  difference 
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W  VELOCITY (CM/S) 


Figure  28.  Time  Series  of  W  Velocity  Component  at  Map  Point  (160,92)  and 


Depth  of  90  Meters  for  Robust  Diagnostic 

between  the  robust  diagnostic  and  the  density  diagnostic  w  component  velocities  is 
discussed  further  below. 
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Figure  29.  Time  Series  of  U  and  V  Velocity  Components  at  Map  Point  (160,92) 
and  Depth  of  65  Meters  for  Robust  Diagnostic. 

The  a,  field  is  allowed  to  change  over  time  in  the  robust  diagnostic.  The  time 
series  in  Figures  24  and  27  show  this  change  for  two  points  in  the  model  domain  Tlie 
initial  a,  field  at  65  meters  depth  for  the  robust  diagnostic  case  is  the  same  as  that  for  the 
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density  diagnostic  (Figure  16).  The  corresponding  field  at  72  hours  is  shown  in  Figure 
30.  The  decrease  in  a,  with  time  shown  at  the  two  points  in  the  time  series  (Figures  24 
and  26)  is  not  uniform  throughout  the  model  domain.  In  general,  areas  of  downweUing 
coincide  with  areas  where  decreases  between  0  and  72  hours,  and  areas  of  upwelling 
coincide  with  areas  where  a*  increases  between  0  and  72  hours  (Figures  31  and  35). 
DownweUing  brings  warmer,  less  dense  water  to  the  65  meter  level  and  upwelling  carries 
colder,  denser  water  to  that  level. 

In  addition  to  the  larger  scale  changes  in  Og  due  to  vertical  advection,  there  has  also 
been  a  tendency  for  adjustment  on  the  smaller  space  scales  to  produce  a  smoother  <Te  field 
at  72  hours.  Some  of  the  smoothing  effect  is  accomplished  by  the  biharmonic  eddy 
diffusion  and  some  by  means  of  the  dynamic  adjustment.  A  scale  analysis  of  the 
thermodynamic  equation  in  the  model  indicates  that  the  parameter  r, 

r  =  {17} 

/  ^ 

where  A  =  biharmonic  difiusion  coefficient,  gives  the  ratio  of  the  two  effects  (i.e.  the 
ratio  of  the  biharmonic  diffusion  term  to  the  change  of  Og  on  the  time  scale  of  adjustment 
(f ')).  Using  the  value  of  A  from  Table  2,  we  find  that  r  =  5  for  L  (wavelength)  =  10  km, 
whUe  r  =  .3  for  L  =  20  km  and  r  =  .06  for  L  =  30  km.  Thus,  the  smoothing  effect  is 
primarily  accomplished  by  the  dynamic  adjustment  for  disturbance  wavelengths  equal  to 
or  greater  than  about  20  km.  Biharmonic  eddy  diffusion  dominates  for  wavelengths  of 
10  km  or  smaUer. 
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The  initial  u  and  v  con^nent  velocity  fields,  and  the  resultant  velocity  vectors  at 
65  meters  depth  for  the  robust  diagnostic,  are  the  same  as  those  for  the  density  diagnostic 
(Figures  17.  18  and  19).  The  corresponding  fields  at  72  hours  are  shown  in  Figures  32, 
33,  and  34.  Comparisons  between  the  u  and  v  velocity  component  fields  at  72  hours  for 
the  density  diagnostic  and  the  robust  diagnostic  show  that  the  robust  diagnostic  fields  are 
smoother  (Figures  20,  21,  32  and  33).  The  smoothness  of  the  u  and  v  component  velocity 
fields  for  the  robust  diagnostic  at  72  hours  is  consistent  with  the  relative  smoothness  of 
the  c,  field  for  the  robust  diagnostic  at  72  hours  as  compared  to  the  initial  Ce  field.  If 
the  noise  in  the  72  hour  u  and  v  component  velocity  fields  in  the  density  diagnostic  case 
is  caused  by  the  effects  of  small  scale  noise  present  in  the  initial  objectively  analyzed 
fields,  the  robust  diagnostic  adjustment  eliminates  these  effects.  Also,  note  that  the 
strongest  u  and  v  component  velocities  at  72  hours  for  the  density  diagnostic  are 
considerably  stronger  than  those  for  the  robust  diagnostic  (Figures  20,  21,  32,  and  33). 
The  resultant  vector  fields  at  65  meters  depth  shown  in  Figures  22  and  34  reflect  this 
difference.  The  density  diagnostic  vector  field  has  a  maximum  velocity  over  100  cm/s 
(Figure  22).  The  robust  diagnostic  vector  field,  on  the  other  hand,  has  a  maximum 
velocity  of  around  75  cm/s,  which  is  closer  to  the  maximum  velocities  stated  above 
observed  by  Huyer  et  al.  (1990)  at  50  meters  depth  during  CTZRl  (Figure  34). 

Conqjaring  the  w  component  velocity  field  for  the  robust  diagnostic  solution  at  72 
hours  (the  initial  w  is  zero  by  definition  for  the  robust  diagnostic  as  well  as  the  density 
diagnostic)  with  the  corresponding  field  from  the  density  diagnostic  solution  shows 
dramatic  difference.  As  noted  in  corutection  with  Figure  28,  the  maximum  upwelling  and 
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Figure  30.  Sigma-9  Field  at  Time  72  Hours  and  Depth  of  65  Meters  for  Robust 


Diagnostic:  Contour  interval  .2  kg/m’. 

downwelling  velocities  for  the  robust  diagnostic  are  an  order  of  magnitude  smaller  than 
those  for  the  density  diagnostic  (Figures  23  and  35).  For  the  robust  diagnostic,  the  time 
series  of  w  component  velocity  at  the  two  points  (Figures  25  and  28)  indicate  that 
downwelling  is  occurring  at  both  points  at  72  hours  at  the  rate  of  approximately  .02  cm/s 
(17  m/d).  For  the  density  diagnostic,  on  the  other  hand,  the  time  series  of  the  w 
component  velocity  (Figures  12  and  14)  indicate  that  upwelling  is  occuring  at  both  point.s 
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CRUISE  CTZR1,  06-12  JULY  1988 


SIGMA  0  HOUR  72  -  HOUR  0.  DEPTH  65  M 


Figure  31.  Sigma-9  Difference  Field  (Time  72  Hours  -  Time  0)  at  Depth  65  Meters 
for  Robust  Diagnostic:  Contour  interval  .2  kg/m’.  Dashed  lines  are 
negative  values. 


at  72  hours  at  the  rate  of  200  to  700  m/d. 

As  indicated  above,  these  vertical  velocities  for  the  density  diagnostic  were  much 
higher  than  the  values  observed  by  Washburn  et  al.  (1990)  and  Kadko  et  al.  (1990).  The 
robust  diagnostic  vertical  velocities  as  represented  by  the  time  series  at  the  two  pouits  luul 
the  w  field  at  90  meters  depth  at  72  hours  (Figures  25,  28,  and  35),  on  the  other  hiuid. 
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Cruise  CTZR1,  06-12  July  1988 


U  hour  72,  depth  65  m 


Figure  32.  U  Component  Velocity  Field  at  Time  72  Hours  and  Depth  of  65  Meters 
for  Robust  Diagnostic:  Contour  interval  20  cm/s.  Dashed  lines  are 
negative  values. 

are  of  the  order  of  20  meters  per  day,  which  are  similar  to  the  above  stated  observed 
values  of  6  to  25  meters  per  day. 

The  advantages  of  the  robust  diagnostic  as  compared  to  the  density  diagnostic  are 
now  clear.  The  robust  diagnostic  solution  comes  to  a  new  balance  much  faster  than  doe.<! 
tlie  density  diagnostic  solution.  In  addition,  at  this  new  state  of  balance  the  robu.st 
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Figure  33.  V  Component  Velocity  Field  at  Time  72  Hours  and  Depth  of  65  Meters 
for  Robust  Diagnostic:  Contour  interval  20  cm/s.  Dashed  lines  are 
negative  values. 


diagnostic  velocity  fields,  horizontal  and  vertical,  arc  much  closer  to  the  observed  values 
than  the  density  diagnostic  Helds.  For  this  reason,  the  balance  of  the  discussion  focuses 
on  the  robust  diagnostic  version  of  the  model. 

Two  minor  drawbacks  of  the  robust  diagnostic  require  further  investigation.  Tl«e 
first,  as  mentioned  above,  is  the  long  oscillations  which  occur  at  the  longer  integration 
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Figure  34.  Resultant  Horizontal  Velocity  Vectors  at  Time  72  Hours  and  Depth  of 


65  Meters  for  Robust  Diagnostic:  Vectors  are  plotted  at  every  4th 


gridpoint. 


tinies  of  the  model.  These  oscillations  could  be  damped  by  shortening  the  damping  time 
scale;  thus,  reducing  the  likelihood  of  tlie  long  wave  synoptic  features  adjusting  with  time. 
This  shortening  of  the  damping  time  scale  could,  however,  reduce  some  of  the  desirable 
geostrophic  adjustment.  Sensitivity  studies,  therefore,  are  needed  to  determine  the 
optimum  damping  time  scale. 
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Figure  35.  W  Component  Velocity  at  Time  72  Hours  and  Depth  of  90  Meters  for 
the  Robust  Diagnostic:  Contour  interval  .02  cm/s.  Dashed  lines 
are  negative  values. 

The  second  drawback  is  the  boundary  noise  seen  in  the  w  component  velocity  field 
at  72  hours  (Figure  35).  The  noise  is  due  to  the  convergence  at  the  boundaries  of  inertia- 
gravity  wave  energy  that  has  rapidly  dispersed  away  from  the  interior  of  the  domain 
during  the  geostrophic  adjustment  process.  This  noise  is  not  as  noticeable  in  the 
horizontal  velocity  fields  because  it  is  primarily  confined  to  the  (smaller)  divergent  part 


55 


of  the  currents.  The  most  effective  method  of  reducing  the  boundary  noise  would  be  to 
increase  the  strength  of  the  sponge,  as  described  above,  or  to  introduce  a  divergence 
damper  in  the  (u,v)  equations  at  the  boundaries.  These  drawbacks,  however,  do  not 
hamper  the  results  of  the  interior  solution  as  presented  here. 


D.  INTERPRETATION  OF  THE  VERTICAL  VELOCITY  HELD  IN  THE 

ROBUST  DIAGNOSTIC  SOLUTION 

One  of  the  most  significant  results  of  the  model  generated  currents  is  the  predicted 
field  of  vertical  velocity.  The  vertical  velocity  is  important  in  determining  the  fate  of 
heat,  salt  and  biological  material  that  is  advected  offshore  in  squirts  and  jets  in  the  coastal 
region.  Comparing  the  w  field  in  Figure  35  with  the  current  field  in  Figure  34  indicates 
that  the  offshore  jet  is  associated  with  sinking  motion,  while  the  onshore  part  of  the  jet 
is  associated  with  rising  motion.  In  this  section  an  attempt  is  made  to  understand  this 
result  in  terms  of  a  quasi-geostrophic  vorticity  budget. 

A  simplified  form  of  the  vorticity  equation  is 


$=-(C 

dt 


{18} 


where  is  the  total  derivative  of  the  relative  vorticity  and  +  V  is  the 
at 

divergence  term  (Holton,  1979;  Pond  and  Pickard,  1983).  The  variation  of  the  coriolis 
parameter  has  been  neglected  for  this  scale  of  motion. 

Amone  et  al.  (1990)  and  Bower  (1991)  describe  the  divergence/upwelling  and 
convergence/downwelling  phenomena  as  observed  in  the  Gulf  Stream  system.  The 
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schematic  shown  in  Figure  36  briefly  summarizes  the  interpretation  by  these  authors.  The 
above  equation  aids  in  analyzing  the  schematic. 

As  seen  from  equation  1 18),  if  divergence  (V-  t^>0)  is  present,  then  the  vorticity 
tendency  of  the  fluid  particle  is  negative,  implying  a  tendency  for  anticyclonic  motion. 
If  convergence  (V-  t^cO)  is  present,  then  the  vorticity  tendency  is  positive,  implying  a 
tendency  for  cyclonic  motion  (Figure  36).  Conversely,  if  a  fluid  particle  flows  from  a 
region  of  negative  vorticity  (a  ridge)  into  a  region  of  positive  vorticity  (a  trough), 
experiencing  a  positive  vorticity  tendency,  then,  according  to  equation  { 18),  it  must  flow 
through  a  region  of  horizontal  convergence. 

From  continuity  and  equation  ( 18),  the  following  relationship  is  true 

K  {19} 

dt  dz’ 


From  equation  1 19)  with  w  =  0  at  z  =  0,  a  negative  vorticity  tendency  implies  upweUing 

(.^<0) ,  and  a  positive  vorticity  tendency  implies  downwelling  (.^>0)  (Pond  and 
dz  dz 

Pickard,  1983).  In  Figure  36,  the  vertical  upweUing  centers  are  shown  upstream  of  a  crest 
in  an  area  of  divergence  and  downweUing  centers  are  shown  downstream  of  a  crest  in  an 


area  of  convergence. 

The  resulting  vertical  velocity  in  the  robust  diagnostic  model  (Figures  35  and  37) 
agree  closely  with  the  panems  of  Amone  et  al.  (1990)  shown  in  Figure  36.  At  72  hours, 
areas  of  downweUing  are  located  downstream  of  a  ridge  in  the  offshore  jet  in  the  'iouth 
central  section  of  the  model  domain  at  65  meters  depth.  The  main  area  of  upweUing  is 
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Figure  36.  Convergence  and  Divergence  Patterns  in  a  Typical  Meander:  Vertical 
upwelling  centers  are  the  dots,  and  downwelling  centers  are  the 
crosses(from  Amone  et  al,  1990). 

located  upstream  of  the  ridge  where  the  flow  in  the  jet  is  directed  onshore  (Figure  37). 

The  vorticity  field  at  the  same  tune  and  depth  (Figure  38)  supports  the  above 
interpretation.  The  strongest  sinking  motion  in  Figure  35  occurs  in  the  region  of  largest 
vorticity  tendency.  This  is  where  fluid  particles  arc  flowing  across  the  vorticity  "front" 
south  of  the  offshore  jet. 

The  initial  pressure  and  vorticity  fields  are  shown  to  illustrate  the  evolution  of  these 
fields  with  time  (Figures  39  and  40).  Note  the  smoothing  which  occurs  in  the  vorticity 
field,  in  particular,  between  the  initial  and  72  hour  times  (Figures  38  and  40).  This 
smoothing  is  the  result  of  the  gcostrophic  adjustments  caused  by  the  robust  diagnostic  as 
discussed  above.  However,  the  basic  pattern  of  positive  vorticity  to  the  left  of  the  Jet 
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(looking  downstream)  and  negative  vorticity  to  the  right  of  the  jet  is  similar  to  the 
vorticity  pattern  of  72  hours. 

Equation  1 19)  can  be  used  to  estimate  the  magnitude  of  upwelling  and  downwelling 
expected  in  the  CTZ  through  scale  analysis.  These  arguments  also  were  used  by  De  Jesus 
(1990).  The  right  hand  side  of  equation  1 19)  can  be  scale  analyzed  as 


dz  H 


where  w.„  is  the  vertical  velocity  at  depth  H,  and  H  is  the  depth  of  interest  (65  meters). 
The  left  hand  side  of  equation  (19)  becomes 


dC 

dt  O 


{21} 


when  scale  analyzed,  where  U  and  L  are  the  velocity  scale  and  length  scale,  respectively. 


Equating  the  right  hand  and  left  hand  sides,  and  solving  for  w  gives 


{22} 

"  L 

where  is  the  Rossby  number. 

Using  the  scaling  estimates  firom  the  advective  and  geostrophic  time  scale  arguments 
described  above 


U  -  50  cm  sec'^ 
L  ~  50  km 
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Figure  37.  Pressure  Field  at  Time  72  Hours  and  Depth  of  65  Meters  for  Robust 


Diagnostic:  Vertical  velocity  pattern  from  Figure  35  shown  with, 
+  =  upwelling,  -  =  downwelling. 


and 


/  -  10  sec 

Tlierefore, 


R  ~  0.1 

O 
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Cruise  CTZR1.  06-12  July  1988 

VORT  hour  72.  depth  65  m 
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o 
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16  32  48  64  60  08  112  128  144  160  176  192  208  224  240 

Distance  (km) 


Figure  38.  Vorlicity  Field  at  Time  72  Hours  and  Depth  of  65  Meters  for  Robust 
Diagnostic:  Contour  interval  1  X  lO"’  sec  '.  Dashed  lines  are 
negative  values. 


w  -  .01  cm  sec  * 


or  about  6  m  day  ',  which  is  similar  to  the  lower  observed  value  of  Washburn  et  nl  ( 1990) 
and  Kadko  et  al.  (1990)  and  similar  to  the  result  shown  in  Figure  35.  These  scaling 


Cruise  CTZR1,  06—12  July  1988 

P  hour  0,  depth  65  m 


Distance  (km) 


Figure  39.  Pressure  Field  at  Time  0  and  Depth  of  65  Meters  for  Robust 
Diagnostic:  Contour  interval  2.5  cm  dynamic  height.  Dashed  lines 
are  negative  values. 


arguments  indicate  that  the  observed  values  and  the  robust  diagnostic  model  values  are 
probably  characteristic  vertical  velocity  values  for  die  CTZ  region.  The  interpretaion  of 
tlie  w  field,  as  that  due  to  the  convergence/divergcnce  pattern  that  changes  the  vorticity 
of  fluid  particles  as  they  move  through  a  meandering  jet  with  alternating  regions  of 
positive  and  negative  vorticity,  is  probably  also  characteristic  of  the  CTZ  region  as  it  is 
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Cruise  CTZR1.  06-12  July  1988 


VORT  hour  0,  depth  65  m 


0  18  32  48  84  80  98  112  128  144  160  176  192  208  224  240 

Distance  (km) 


Figure  40.  Vorticity  Field  at  Time  0  and  Depth  of  65  Meters  for  Robust 
Diagnostic:  Contour  interval  1  X  10’’  sec'*.  Dashed  lines  are 
negative  values. 

of  the  Gulf  Stream  (Bower,  1991). 


( 


E.  ADCP  COMPARISONS 

Tlie  resultant  velocity  field  for  the  robust  diagnostic  at  time  72  hours  and  65  meters 
depth  shows  a  relatively  smooth  southward  to  southwestward  flowuig  jet  tiuough  tlie 
center  of  the  model  domain  (Figure  34).  The  magnitude  of  the  velocity  following  the  jet 


63 


is  fairly  steady,  and  it’s  direction  changes  smoothly.  The  ADCP  velocity  vectors  for 
CTZRl  at  50  meters  depth  (Figure  41),  on  the  other  hand,  indicate  the  instantaneous 
observed  velocity  field  is  not  quite  so  smooth  as  indicated  in  the  model.  The  general 
shsqTe  of  the  jet  in  the  observed  field  is  reflected  in  the  model  field.  The  magnitude  of 
the  velocities  in  the  observed  field,  however,  change  spatially  to  a  much  greater  degree 
and  with  more  irregularity  than  in  the  model.  This  is  probably  to  be  expected  since  the 
ADCP  data  represent  instantaneous  currents,  while  the  model  currents  represent 
geostrophically  balanced  currents  estimated  from  a  fairiy  coarse  sample  grid.  Also,  the 
maximum  observed  velocities  are  in  the  southwestern  comer  of  the  cmise  domain  where 
the  jet  is  farthest  offshore.  In  the  model,  maximum  velocities  for  the  jet  are  indicated  in 
the  northeastern  part  of  the  domain,  closer  to  the  shore. 

The  smoothing  which  occurs  because  of  the  robust  diagnostic  seems  to  eliminate 
some  of  the  [^ysical  small  scale  disturbances  in  the  velocity  field  as  well  as  those 
introduced  by  the  objective  analysis  as  discussed  above.  Vertical  cross-sections  of  the 
velocity  field  (below)  also  indicate  diat  dre  model  does  not  represent  well  the  smaller 
scale  characteristics  of  the  observed  velocity  field.  This  is  to  be  expected,  to  some  extent, 
because  the  model  currents  represent  a  balanced  state,  while  the  ADCP  data  are 
instantaneous  currents. 

Veitical  cross-sections  of  normal  velocity  from  the  model  domain  are  compared  to 
the  corresponding  ADCP  normal  velocity  cross-sections  from  the  cruise  (Figure  42).  The 
cross-sections  corresponding  to  cruise  leg  B  in  the  northeastern  comer  of  domain  sliovv 
that  the  observed  jet  (Figure  43)  is  smaller  in  venical  and  horizontal  extent,  and  has  a 
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leg  C  (Figure  45)  corresponds  closely  with  that  developed  by  the  model  (Figure  46),  but 
the  observed  intensity  of  the  jet  is  not  as  great  as  that  developed  by  the  model.  The 
cross-sections  for  cruise  leg  D  have  similar  characteristics  to  those  of  cruise  leg  C(Figures 
47  and  48).  Cruise  legs  E,  F  and  G  cross-sections,  however,  show  an  increasing 
agreement  in  the  magnitude  as  well  as  the  position  of  the  observed  jet  and  the  modeled 
jet  (Figures  49,  50,  51,  52,  53  and  54).  In  general,  note  that  the  model  velocity  cross- 
sections  for  all  legs  are  much  smoother  than  those  for  the  ADCP  velocities.  Smaller  scale 
disturbances  do  not  appear  in  the  model  cross-sections  because  of  the  balanced  nature  of 
the  currents  that  are  represented,  and  the  lower  resolution  of  CTD  stations  used  to 
initialize  the  model.  The  discrepancy  between  model  and  ADCP  data  below  the  surface 
jet  (i.e.  below  200  m),  however,  cannot  be  attributed  to  smaller  scale  imbalanced  flows 
in  ADCP  data  because  the  currents  at  that  depth  should  be  in  dynamic  balance. 

There  are  several  possible  reasons  why  the  ADCP  and  model  currents  are  in  closer 
agreement  farther  offshore.  Topogrqthic  effects,  which  are  neglected  in  the  model,  could 
produce  the  increasing  agreement  between  observed  velocities  and  modeled  velocities 
with  progression  from  northeast  to  southwest  through  the  model.  Boundary  effects  in  the 
model  could  also  play  a  part,  especially  in  the  northeast  comer  of  the  model  domain 
where  inflow  of  the  jet  is  occurring.  In  addition,  the  boundary  probably  affects  the  model 
cross-section  corresponding  to  cruise  leg  G  (Figure  54).  Note  the  reverse  flow  which 
occurs  in  the  extreme  southeastern  part  of  this  leg  which  docs  not  occur  in  the  observed 
velocity  field  (Figure  53).  More  rigorous  statistical  studies,  and  robust  diagnostic  studies 
with  topography,  wind  and  heat  flux  forcing,  and  barotropic  currents  included,  are 
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Figure  42.  Model  Domain  Cross-Sections,  Cruise  CTZRl:  Solid  lines  indicate 
model  cross-sections.  Diagonal  rows  of  dots  are  cruise  legs  A 
through  G  from  northeast  to  southwest. 

required  before  a  better  interpretation  of  the  comparison  between  ADCP  data  and  model 


currents  can  be  made. 


CTZ88/GRID-IB 

ADCP  VELOCITY  (CM/S) 

10  KM 

7/7:1223 

7/8:430 

Figure  43.  ADCP  VelocKy  Cross-section  Cruise  CTZRl  Leg  B:  Arrows  are 
boundaries  of  corresponding  model  cross-section.  Dashed  lines 
indicate  northeastward  flow.  Northwest  is  to  left.  Contour  interval 
5  cm/s 
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Figure  45.  ADCP  Velocity  Cross-section  Cruise  CTZRl  Leg  C:  Conventions  and 
contour  interval  are  same  as  in  Figure  43. 
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Figure  46.  Model  Velocity  Cross-section  Cruise  CTZRl  Leg  C:  Dashed  lines  are 


northeastward  flow.  Northwest  is  to  left.  Contour  interval  5  cm/s. 
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Figure  47.  ADCP  Velocity  Cross-section  Cruise  CTZRl  Leg  D:  Conventions  and 


contour  interval  are  the  same  as  in  Figure  43. 


12 


Depth  (m) 


Cruise  CTZR1,  06—12  July  1988 
norm  field  at  hour  72  from  (112,156)  to  (188,80) 


NW-SE  Distance  (km) 


Figure  48.  Model  Velocity  Cross-section  Cruise  CTZRl  Leg  D  :  Dashed  lines  are 
northeastward  flow.  Northwest  is  to  left.  Contour  interval  5  cm/s. 
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Figure  49.  ADCP  Velocity  Cross-section  Crube  CTZRl  Leg  E:  Conventions  and 
are  contour  interval  are  the  same  as  in  Figure  43. 


Cruise  CTZR1.  06-12  July  1988 


norm  field  at  hour  72  from  (88,128)  to  (156,60) 


NW-5E  Distance  (bm) 


Figure  50.  Model  Velocity  Cross-section  Cruise  CTZRl  Leg  E:  Dashed  lines  are 
northeastward  flow.  Northwest  is  to  left.  Contour  interval  5  cin/s 
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Cruise  CTZR1,  06-12  July  1988 


norm  field  at  hour  72  from  (44,120)  to  (132,32) 


NW— 5E  Distance  (km) 


Figure  52,  Model  Velocity  Cross-section  Cruise  CTZRl  Leg  F:  Dashed  lines  are 
northeastward  flow.  Northwest  is  to  left.  Contour  interval  5  cin/s. 
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CTZSa/GRID  IG  I - 

ADCP  VELOCITY  (CM/S)  10  LM 


Figure  53.  ADCP  Velocity  Cross-section  Cruise  CTZRl  Leg  G:  Solid  lines  are 
southwestward  flow.  Northwest  is  to  left.  Contour  interval  5  cin/s. 
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Cruise  CTZR1,  06—12  July  198b 

norm  field  ot  hour  72  from  (12,92)  to  (100,4) 


Figure  54.  Model  Velocity  Cross-section  Cruise  CTZRl  Leg  G:  Arrows  are 
boundaries  of  corresponding  ADCP  cross-section.  Dashed  lines  are  noitteastwaid  flow. 
Northwest  is  to  left.  Contour  interval  5  cm/s. 
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rv.  SUMMARY  AND  CONCLUSIONS 

The  "interface  matching"  formulation  of  determining  the  deep  solution  for  the 
disturbance  gives  reasonable  results  when  input  into  the  primitive  equation  model  used 
in  this  study.  No  comparisons,  however,  were  made  between  these  results  and  those 
which  would  be  obtained  using  the  "mean  conserving"  or  the  "energy  conserving" 
formulations.  A  deep  profile  comparison  of  the  "interface  matching",  the  "mean 
conserving",  and  the  "energy  conserving"  formulations  with  observed  data  is  needed  to 
determine  which  method  gives  results  which  most  closely  resemble  the  observed  vertical 
ocean  structure.  In  addition,  a  comparison  is  needed  between  the  results  obtained  by  the 
exact  method  and  those  obtained  by  the  linear  least  squares  method  of  solving  the 
minimization  problem  set  up  between  equations  { 1 )  and  (2)  as  described  in  Chapter  II. 

The  robust  diagnostic  version  of  the  model  gives  better  results  than  the  density 
diagnostic  version  of  the  model  for  this  case.  The  robust  diagnostic  version  of  the  model 
allows  a  mutual  dynamic  adjustment  to  occiv  between  the  Oe  and  velocity  fields. 
Consistent  with  the  theory  of  geostrophic  adjustment,  the  smaller  scale  features  in  the 
field  are  altered  the  most  by  this  adjustment  process.  The  improvement  of  the  robust 
diagnostic  over  the  density  diagnostic  versions,  therefore,  will  probably  hold  true  for  the 
other  CTZ  cases  as  well.  The  accuracy  with  which  the  robust  diagnostic  version  of  the 
model  diagnoses  the  magnitude  and  pattern  of  the  observed  vertical  velocity  associated 
with  jets  in  the  CTZ  has  not  been  accurately  determined.  However,  the  model  results  that 
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indicate  horizontal  convergence  and  sinking  motion  of  the  order  of  5-20  m/day  all  along 
the  offshore  directed  part  of  the  jet  are  certainly  consistent  with  the  few  observed 
estimates  that  are  presently  available.  A  sensitivity  study  to  determine  the  optimum 
damping  time  scale  for  the  robust  diagnostic  for  reducing  advective  effects  is  desirable, 
but,  as  stated  in  Chapter  in,  the  fact  that  advective  effects  were  apparent  in  this  particular 
robust  diagnostic  version  did  not  adversely  affect  the  results. 

Comparisons  of  ADCP  with  robust  diagnostic  model  velocity  vertical  cross-sections 
indicated  that  small-scale  features  of  the  observed  velocity  Held  were  not  reflected  in  the 
model  velocity  field.  This  smoothing  of  small-scale  features  by  the  robust  diagnostic  was 
seen  in  all  fields,  horizontal  as  well  as  vertical.  The  smoothing  eliminated  some 
undesirable  noise  resulting  from  the  multi-layer  objective  analysis  scheme,  but  also 
eliminated  some  physical  features. 

The  vertical  cross-section  comparison  between  ADCP  velocities  and  robust 
diagnostic  model  velocities  also  showed  that  agreement  between  the  two  increased  from 
northeast  to  southwest  in  the  model  domain.  In  addition,  ADCP  velocities  were  higher 
in  the  jet  when  it  was  farthest  offshore  in  the  southwest  comer  of  the  model  domain  while 
model  velocities  were  higher  in  the  jet  where  it  was  closer  to  shore  in  the  northeast 
comer.  A  hypothesis  for  this  difference  is  the  model’s  neglect  of  topographic  effects, 
wind  and  heat  flux  forcing.  An  interpretation  of  robust  diagnostic  model  results  with 
these  effects  added  is  required  to  test  this  hypothesis. 
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